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An axially symmetric solution of Laplace's equation in oblate spheroidal coordinates is 
found, which may be used as the gravitational potential about an oblate planet. This 
potential, which makes the Hamilton-Jacobi equation for a satellite orbit separable, has an 
expansion in zonal harmonics in which the amplitudes of the zeroth and second harmonics 
can be adjusted to agree exactly with the values for any axially symmetric planet and a fourth 
harmonic which then agrees approximately with the latest value for that of the earth. 
The net result is therefore a reduction of the problem of satellite motion to quadratures, with 
use of a potential field that is much closer to the empirically accepted one for the earth than 
any heretofore used as the starting point of a calculation. It may thus be possible to do the 
gravitational theory of a satellite orbit very accurately without use of perturbation theory. 
The method can take into account a first harmonic in the potential, in case observations are 
reduced to a center which does not coincide with the center of mass of the planet. 

1. Introduction 

The purpose of this paper is to develop a method of solution of the central problem of 
satellite astronomy, viz, the theory of the effect of the oblateness of a planet on the orbit of any 
satellite, bound or unbound, but near enough to the planet so that the forces of other bodies 
may be neglected. (A preliminary note has appeared in [I]. 2 ) 

To do so, let us consider the gravitational potential V outside an axially symmetric oblate 
planet of mass M. We call the axis of symmetry the polar axis; this will ordinarily be indistin- 
guishable from the axis of rotation, except for entirely negligible periodic effects. We also 
choose some point on the polar axis as the "center" or origin of a coordinate system and let r 
be the corresponding position vector of a field point. With r = |r|, the declination as the com- 
plement of the angle between r and the poiar axis, G the gravitational constant, M the mass of 
the earth, and \x = GM ) the potential V(r, 6) can be expressed by means of an expansion in 
spherical harmonics : 

F=- M r- 1 [l-g Jn(R/r) n P n (sin »)]• (1) 

Here R is the equatorial radius, i.e., the radius of a section through perpendicular to the polar 
axis. The notation J n for the coefficients was introduced by Merson and King-Hele [2]. 
For points on the polar axis (6=ir/2) at the distance R from the center, the ratio of any two 
terms in (1) is equal to the ratio of the corresponding J's. 

If the planet has symmetry with respect to the equatorial plane, all the odd harmonics 
drop out; i.e., Ji=J 3 =J 5 =. . . = 0. (In any case, even if no such symmetry exists, the first 
harmonic drops out, i.e., J x vanishes, if is chosen at the center of mass.) The coefficients 
J 2 and J A are given by 

J 2 =2J/>3=(I p -I t )/MR 2 , (2.1) 

J 4 =-8Z>/35, (2.2) 

where the quantities J and D were introduced by Jeffreys [3], and where I v and I t are the 

1 This work was supported by the U.S. Air Force, through the Office of Scientific Research of the Air Research and Development Command. 

2 Figures in brackets indicate the literature references at the end of this paper. 
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moments of inertia of the planet about polar and transverse axes through 0. 
with taken at the center of mass, Jeffreys [3] gives the values 



For the earth, 



J=(1637.5)10- 



Z>=(10.7)10- 6 , 



while King-Hele and Merson [4] derive from recent satellite observations the values 
J=(1624.6±0.3)10- 6 , £>=(6±1)10- 6 , J 6 =(0.1±1.5)10" 6 . 



(3.1) 



(3.2) 



It should be noted that until very recently it has always been assumed that the earth has 
symmetry with respect to its equatorial plane — strictly speaking, that the geoid, or smoothed 
earth, has such symmetry. That is, all the odd J's have been assumed to vanish. Recently, 
however, there has appeared some evidence [5] for an observable asymmetry with respect to 
the equatorial plane. Such evidence, if confirmed, would have the following results. If 
is taken at the center of mass C, then J x would still vanish, but J 3 , J 5 , . . . would have small 
nonvanishing values. If, on the other hand, is taken, either by error or by design, as non- 
coincident with C, but still on the polar axis, then we should have Jit^O, with J 3 , J 5 , . . . still 
small but nonvanishing. The only conceivable difficulty that might arise with this latter 
choice for the center has to do with the inertial forces acting on the satellite. This difficulty 
is discussed and removed in the next section. A reason for deliberately choosing as non- 
coincident with C might be to represent any equatorial asymmetry largely by means of a first 
harmonic, and such a procedure might work if the third and higher odd harmonics then turn 
out to be negligible in comparison with the first. 

In the present paper we first derive a solution V L of Laplace's equation, in oblate spheroidal 
coordinates, having axial symmetry and leading to exact separability of the Hamilton-Jacobi 
equation for the motion of a satellite. It contains three adjustable constants, bo, bi, and c, 
which may be so chosen that the expansion of V L in spherical harmonics is 

V L =-^ 1 [l^J 2 (B/r) 2 P 2 (sm$)+Ji(E/ryP,(sm . . 

^J 1 (RJr)P 1 (sin 6) + J 1 J 2 (R/ryP 3 (sm 0)- JiJfOZ/r) 5 ^ (sin $) + . . .], (4) 

where J 2 is given by (2.1) and where J x is the same as the J x in (1), viz, 

J 1 =-d/R. (4.1) 

d being the northerly displacement of the center of mass C from 0. 

Comparison of (4) with (1) shows that the fit is exact for the zeroth harmonic, — /xr _1 , 
and for the first and second harmonics. It would be exact through all harmonics if the 
planetary distribution of mass happened to be such that 



e/4 — t/ 2 , J 6 — <S2) 






(5.1) 
(5.2) 



To obtain some idea of the adequacy of the potential (4) let us consider only the even 
harmonics. From (3) and (4) we may construct the following table of their amplitudes, with 
signs. 





Table 1 






Source 




Harmonic 






(Zeroth) 106 


-106J" 2 


-IO6J4 


-106J6 


Jeffreys [3]. 


10 6 

106 
106 

106 


-1090 

-1090 
-1080 

-1080 


2.45 

1.19 

1. 37±0. 23 

1.17 




From fit of 60 and c to Jef- 
freys' J __. __ __ 


-(1.30)10-3 


King-Hele and Merson [4].. 

From fit of bo and c to their 

J 


-0. l=bl. 5 

- (1. 26) 10~ 3 
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The fit of the zeroth and second harmonics is necessarily exact. If we fit to Jeffreys' J", 
the resulting fourth harmonic has an amplitude amounting to 49 percent of that given by 
Jeffreys; the resulting error of (4) then amounts to 1.26 ppm, corresponding to an error of 1.26 
parts in 1100 for the oblateness potential. If we fit to King-Hele and Merson'se/, the resulting 
fourth harmonic has an amplitude amounting to (86 ±16) percent of their value for the fourth 
harmonic, so that the agreement is within the observational error of their determination. 
Even the sixth harmonic as fitted also lies within their estimates, which are, however, \«r\ 
hazy. In any event the resulting error of (4) then amounts to only 0.20 ppm, corresponding 
to an error of only one part in 5500 for the oblateness correction 

The potential function that we find in this paper, leading to Hamilton-Jaeobi separability, 
is thus much closer to the actual potential, at every point of space, than any heretofore used 
as the starting point of a calculation. If the results of King-Hole and Merson [4] are con- 
firmed, it may be possible to do orbit theory for unretarded satellites very accurately without 
perturbations of the solution of the Hamilton-Jaeobi equation. If not, there is still a good 
chance of minimizing the use of perturbation theory, by distributing over the zeroth and 
second harmonics the secular effect of the fraction of the fourth harmonic t hat is not accounted 
for. 3 

2. Reference System 

As a reference system we first choose a right-handed system of rectangular coordinates, 
with origin at the center 0, with Oz pointing towards the planet's north celestial pole and with 
Ox pointing towards the planet's vernal equinox. The directions of the axes are then fixed in 
an inertial system. For the computation of satellite orbits a slow precession of the axis of 
symmetry about the celestial polar axis may be disregarded, if the precession is sufficiently 
slow and if the angle of the cone of precession amounts only to about 0.1 ", as in the case of 
the earth [3, p. 202]. Then Oz may be considered coincident with the axis of symmetry. 

The forces per unit mass acting on a satellite are then, relative to such a system and in 
the absence of drag : 

Fi= gravitational field of the planet=— VV, (6.1) 

F 2 =vector sum of the gravitational fields of the sun and of any moons that the planet 
may have, plus the inertial forces on the satellite produced by acceleration of the center of 
mass of the planet by the sun and the moons, (6.2) 

F 3 = inertial force on the satellite produced by acceleration of the point in the planet, 
relative to its center of mass C, = — 6, where 8 is the vector from C to 0. (6.3) 

If in (1) we disregard all higher harmonics except the second and calculate — W, we find 
that the ratio of the amplitude of the oblateness force F to the total gravitational force is of 
the order 4.5J 2 ~ 1/200. Thus, for the earth 

iF |<5 cm/sec" 2 . (7.1) 

Furthermore [6], 

|F 2 |«10- 4 cm/sec- 2 (7.2) 

and 

|6| = (2tt/T 2 ,) 2 5 sin a« (1.4)10- 20 S sec" 2 (7.3) 

for a period T p of precession «433 days [3, p. 202] and for an angle a of the cone of precession 
«0.1". Thus 

F Z /F 2 ~(1A)10- U 8 (<5 in kilometers), (7.4) 

so that F 3 is utterly negligible compared to F 2 , for any conceivable displacement 8 of from C. 
In turn, F 2 amounts to less than 1/50,000 of the maximum oblateness force. 



The author is indebted to Dr. Boris Garfinkel for this suggestion. 
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We are thus left with Fi as the only appreciable force in the absence of drag, and we may 
choose our coordinate origin either at the center of mass C or at any near point on the polar 
axis that leaves J 2 essentially unchanged but simplifies the odd harmonic part of (1). 

3. Oblate Spheroidal Coordinates 

We next introduce the oblate spheroidal coordinates £, rj, and <£, defined by 

x=c[(?+l) (1-?7 2 )P cos (0<K ») , (8.1) 

y=c[(e+l)(l-r } 2 )]*sm$ (-15*^1), (8.2) 

s=<*l, (8.3) 

where c is a distance to be fixed later. (For many of the following calculations, the book by 
Morse and Feshbach [7] is a very useful reference.) The surfaces £= constant are oblate 
spheroids; sections perpendicular to the £-axis are circles and meridian sections are confocal 
ellipses of semiaxes c(£ 2 +l)^> and c£, the focal separation being 2c. The locus of the foci is 
thus a circle of radius c in the plane 2=0. The surfaces 7?=constant are one-sheet hyperboloids 
of revolution, asymptotic to cones with vertices at the origin, and the surfaces with right 
ascension <£=constant are planes containing the polar axis. Note that rj= ± 1 along Oz or — Oz 
respectively and that for sufficiently large values of r= (x 2 -\-y 2j rZ 2 )^ the coordinates £ and t\ 
behave like £~r/c and Tj^sin 6, where 6 is the declination [7, p. 1292]. The plane z=0 consists 
of the regions 77=0 and £=0, outside and inside the focal circle, respectively. 
The element of arc ds is given by 



where [7, p. 1292] 



ds 2 =hW+h 2 dy 2 +hW 2 , (9) 

h^cie+^He+iy* (10.1) 

h 2 =c(e+n 2 )Hi-n 2 P (10.2) 

hs=c({t 2 +mi-y 2 )>- (10.3) 
4. Hamilton-Jacobi Equation 

Henceforth all quantities proportional to the mass of the satellite are taken per unit mass 
of the latter. By (9) the kinetic energy of the satellite is then given by 



so that 



T=\ (hlHhW+hU 2 ), (11) 



t=\ (hrYi+KYv+Wpl), (12) 



where the generalized momenta are 



S(%, v, <t>) being the action function. 



n=hU = dS/dt (13.1) 

p,=hlv=Z>Sldr, (13.2) 

p^=h i 3 4>=dS/d4>, (13.3) 
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For the case of axial symmetry V=V(£, rj) and the Hamiltonian H becomes 

H = (2c2)-1 [|+? *+£$ *+ v+m-*)} *- v <*>* (14) 

The Hamilton-Jacobi equation is then 



(f 2 + 1) (5^/Sf ) 2 + (1 -H 2 ) (d^/diy) 2 + (ga^^iL^gj ^SI^r + 2c^e+^)V{^)=2c\e + n 2 )ccu 



where a x is the constant energy. 

To try to separate variables in (15), we express S as 



(15) 



S=S 1 (i-)+S 2 (v)+S 3 (4>), (16) 

whereupon (15) becomes 

(? 3 +l)S; 2 +(l-, 2 )^ 2 +(^-^ I )^ 2 +2c 2 (J 2 +, 2 )F(f,,)=2c 8 (| 2 + 1? 2 )a 1 . (17) 

Here 

S'^dSJdt, S' 2 =dS/ 2 d v S' 3 =dS z /d<f>. (18) 

From (17) it follows that S' 3 2 is expressible as a function of £ and 77 alone, so that 

S3=a 3 =constant. (19) 

Further inspection of (17) shows that the variables £ and 77 can be separated if and only 
if the potential F(£, rj) has the form 

V&r } ) = (¥+r) 2 )-V(Q+<jM]. (20) 

The question then arises: Is it possible to find a solution of Laplace's equation, of the form (20), 
thai will adequately represent the known facts about the earth's potential? As indicated in 
the introduction, the answer is essentially yes. We accordingly devote the next section to 
finding all the possible solutions of Laplace's equation that have this form. 

5. Solutions of Laplace's Equation That Lead to Hamilton-Jacobi Separability 

For axial symmetry, Laplace's equation in oblate spheroidal coordinates becomes [7, p. 
1292] 

& < *- i $KI>-<F <21) 

Insertion of (20) into (21) then shows that the functions /(£) and g(rj) must satisfy the equation 

f 2 (« 2 +l)/"(f)-2f(| 2 +2)y'(J)+2(| 2 +2)/(|)+, 2 (l-, 2 ) S r"(,)+2,(, 2 -2) 3 '(,)-2(, 2 -2) 5 (,) 

+M(£,„)=0, (22) 
where 

M(S, i)mifF(&-?0h) (23) 

^'(?) = (£ 2 +l)/"(!)+2£/'(£)-2/(£) (24.1) 

G( v ) = ( v °-l)g"(T,)+2 v g'(t,)-2g(r,). (24.2) 

By (22) a solution of (21) of the form (20) can exist only if M(|, rj) is separable. Now the neces- 
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sary and sufficient conditions that M(£, 17) be separable are that F(£) and #(77) have the forms 

^(5)=4B+4(7£ 2 (25.1) 

G f (i 7 ) = -44+4CV, (25.2) 

where -4, 5, and C are constants and where the signs and the factor 4 are chosen with a view to 
later convenience. The conditions are obviously sufficient; to show that they are also necessary, 
note that if M(£, 77) is separable, there must exist functions <£(£) and \p(rj) such that 

i?F(i)-eG(v)=*(*)-+i&. (26) 

On first placing rj=0 in (26) and then £=0, we find respectively 

<t>& =+(0) + G(0)e (27.1) 

^) = <K0) + ^W. (27.2) 

Then, on inserting the expressions (27) into (26), we find 

V 2 [Ftt)-Fm-?[G(v)-G(0)] = <l>(0)-H0). (28) 

Placing either £=0 or 77=0 in (28) now shows that <£(0) = ^(0), 
so that 

r[FU)-F(0)]=riG(v)-G(0)] = ^C, (29) 

where C is a constant, by a familiar argument. Thus 

F(0 = F(0)+4:Ce (30.1) 

G(ri) = G(0)+4Cr?, (30.2) 

which have the general forms indicated in (25), so that the latter are also necessary. 

Comparison of (24.1) and (25.1) and of (24.2) and (25.2) shows that the functions /(£) and 
g(rj) must satisfy 

(t 2 +l)1"+2tf'-2f-4C?-4tB=0 (31.1) 

(r?-l)g" + 2rig' — 2g-4:Cri 2 + 4A = 0. (31.2) 

By (23) and (25) it now follows that 

M(Z, V ) = 4A?+±B V \ (32) 

Insertion of (32) into (22) then shows that the functions /(£) and g(rj) must also satisfy 

V 2 (v 2 -l)9 f '(v)-2v(v 2 -2)g , (v) + 2( v 2 -2)g(rj)-4Br ] \ (33) 

The usual argument about separability then shows that each side of (33) is equal to a constant 
iK. Thus/(£) and #(r/) must also satisfy 

r(S 2 +l)/"-2!a 2 +2)/'+2($ 2 +2)/+4^f 2 -4K=0 (34.1) 

V 2 (v 2 -l)g"-2r,tf-2)g' + 2(r, 2 -2)g-±B v >—4K=0. (34.2) 
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We may avoid solving second-order differential equations, as follows: To eliminate/'', 
multiply (31.1) by £ 2 and subtract the resulting equation from (34.1). To eliminate g", multiply 
(31.2) by 7j 2 and subtract the resulting equation from (34.2). The resulting first-order equations 
are the linear ones 

lf-j=(e+l)- l [C?+(A+B)t-K\ (35.1) 

r,g'-g=(v 2 -l)- 1 [Cv i -(A+BW-K\. (35.2) 

Their most general solutions are 

M)=K+boZ+Ce+(A+B-C+K)l; tan-^ (36.1) 

g(v) = -K-b lV + Cv 2 +l (A+B-C+K)r, In ^±2 (36.2) 

where b and b x are constants of integration. Insertion of (36.1) into (31.1) now shows that 

B=A. (37) 

Then, with B=A } insertion of (36.2) into (31.2) simply yields a check. 

Thus 

M) = K+b£+Cl 2 +(2A-C+K)l tan- 1 ? (38.1) 

g( v )= -K-b lV +C v 2 +l(2A-C+K)r J In j±J ( 38 .2) 

It is easy to verify that these solutions for/© and g(rj) satisfy both sets of second-order equa- 
tions, viz, (31) and (34), as they should. Moreover, since/© and g(rj) must satisfy (35) and 
since the expressions (36) are the most general solutions of (35), it follows that there are no 
other solutions for/© and g(rj). 

From (38) it now follows that all axially symmetric solutions of Laplace's equation that 
have the form (20) in oblate spheroidal coordinates are expressible as 

V(^)=C+(^+ v 2 )- 1 [b ^-b l v+b i (2^ ttin-^+r, m^)]> (39) 

where we have replaced ^(2A—C J rK) by b 2 . The Hamilton-Jacobi equation is separable, 
when the potential is a solution of Laplace's equation, if and only if the potential has lliis 
form in oblate spheroidal coordinates. 

If as usual we take V to be zero at infinity, the constant C vanishes. We are thus left 
with three explicit adjustable constants, b , b iy and b 2 . It is at once clear, however, that b 2 
must vanish if the solution is to represent the potential outside a planet, for the logarithmic 
term has a singularity everywhere on the polar axis, for which rj= ± 1. 

We thus finally arrive at a solution of Laplace's equation which makes the Hamilton- 
Jacobi equation separable and which contains two explicit adjustable constants, 6 and b u 
and one implicit adjustable constant c, the latter appearing in the equations of transformation 
from oblate spheroidal to rectangular coordinates. This solution is 

V=(¥+V 2 )- 1 (bo!i-b l r } )=b Re ft+i^+ft, Im (S+jrj)- 1 . (40) 

The reader may readily verify that (%+jy)' 1 is a solution of the Laplace equation (21) and thus 
that (40) is a solution. (We here use j instead of i for the imaginary unit, saving i for use in 
later papers as orbital inclination.) 
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6. Expansion of the Potential in Spherical Harmonics 

From (8) we obtain 

^-rj 2 =r 2 /c 2 -l, (41.1) 

fr=sr/c, (41.2) 

where r is the distance from the origin to the satellite and where 

s = sin 6, (41.3) 

6 being the declination of the satellite. Then 

(i;+jv) 2 =¥-r) 2 +2j£ri=r*c- 2 -l+2jsrc- 1 (42.1) 

= r 2 c- 2 (l-2hs + h 2 ), (42.2) 

where 

h = -jcr-\ (42.3) 

Then 

^+jv)~ l = cr-\l-2hs+h 2 )-^. (43) 

Since we shall choose c<<r, it follows that |A|<1 and (1 — 2hs+h 2 )~^ is then the generating 
function for the Legendre polynomials P n (s). Thus 

(H+jvy^cr- 1 ± h n P u (s)=cr-' S (-j) n c n r-»P n (s), (44) 

from which there follow the relations 

(i 2 + V 2 )~ 1 ^^^+jr ) )- 1 = cr- 1 [l-c 2 r- 2 P 2 (sm <9) +c 4 r~ 4 P 4 (sin 0)- . . .] (45.1) 

(¥+V 2 )- 1 V=-I^+h)~ 1 = cr- 1 [cr- 1 P l (sin 0)-c 3 r- 3 P 3 (sin 0) + c 5 r- 5 P 5 (sin $) — .. .]. (45.2) 
From (40) and (45) it then follows that 

V= b cr~ ! [1 -c 2 r~ 2 P 2 (sin 6) +c 1 t^P a (sin 0) — . . .]— ijcr-^cf 1 ?! (sin 0)— eV- 3 P 3 (sin 0) + . . .] 

(46) 

7. Evaluation of the Adjustable Constants 

Comparison of (46) with (1) shows that to fit the zeroth harmonic we must choose 

hc=-^ (47.1) 

that to fit the second harmonic we must then choose 

c 2 =R 2 J 2 , (47. 2) 

and that to fit the first harmonic we must choose 

6 1 c 2 =-MSe/ 1 =^. (47.3) 

That is all the flexibility that (46) permits; the higher harmonics are then all fixed, as indicated 
in the relations (5), which follow from (46). 

These choices lead to the form (4) for the potential, which has been amply discussed in 
the introduction. For i?=6378. 388 km we find from (47. 2) and table 1 : 

c=210. 7 km (Jeffreys), (48. 1) 

c=209. 9 km (King-Hele and Merson). (48. 2) 
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Although eqs (8) appear to place this constant c in so central a position in 1 lie theory as to make 
numerical predictions depend strongly on its accuracy, it is not really so. Actually we shall 
see that a knowledge of the value of c is needed only to find the deviations of the orbit from a 
conic section. 

From (40) and (47) we obtain as our working expression for V, in oblate spheroidal co- 
ordinates : 

v=- f jLC-\e+r 1 *y i ti+vVc)- (49) 

8. Formal Solution of the Hamilton-Jacobi Equation 
Insertion of (19) and (49) into (17) now results in the fully separated equation 

(eM-l)S{^-(?+l)-y-2AU*-^ (50) 

By the usual argument, each side of (50) is equal to a constant k. 

At this point it is necessary to discuss complete and incomplete orbits. We define a 
complete orbit as the orbit that would be traversed by the satellite when the time t is allowed 
to run from — °° to + °°, on the assumption that the equations defining the orbit hold for all values 
of t. In the cases of a planetary satellite or of a meteor coining in from infinity with its orbit 
being bent around the planet without hitting it, the actual physical orbit is complete. Any 
projectile fired from a gun and returning to earth, however, exemplifies an incomplete bound 
orbit. Any missile which escapes from the earth exemplifies an incomplete unbound orbit, as 
does any meteor with energy ai>0 which hits the earth. Note that the complete orbit may or 
may not intersect the planet. If it does, we define it as that complete orbit which would be 
described by the satellite, in the absence of planetary matter, but in a gravitational field cor- 
responding to (49). For r^>c, this field can be represented by the multipole field (46). 

By (50) it follows that 



k+al- 



-(l-^-^U-i? 2 )- 1 ^- 



-2 i u5i7 + 2a 1 cV. 



(51) 



Thus if the plane 77=0 is reached by the complete orbit 

k+a 2 3 ^0. 



(52) 



In the case where 5 is chosen to be zero, the plane 77 = is the ordinary equatorial plane through 
the center of mass of the planet. Thus, if this plane is reached by the complete orbit, (52) 
is satisfied. 

It is clear that (52) is only a necessary condition that the complete orbit reach the plane 
77=0. It is difficult to find a corresponding sufficient condition. By analogy with the Kepler 
problem (r _1 potential), we may suspect that when 5=0 the equatorial plane will always be 
reached by the complete orbit of any bound satellite. If 5^0, however, the plane 77 = may 
not be reached if the orbit is very close to the equator. This latter circumstance is illustrated 
in figure 1. 



I 



POLAR AXIS 



PROJECTION OF ORBIT 
ON PLANE THROUGH 
POLAR AXIS 




PLANE 17 = 
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The difficulties in finding a condition sufficient for 77=0 to be reached are illustrated 
by the existence of cases where the complete orbit may not reach either the region 77=0 
outside the focal disk or the region £=0 inside it. For example, if a particle is dropped from 
a point on the polar axis into a matter-free region satisfying (49), so that a 3 =0, 770= ±1, Vo=0 } 
it follows (with 5 = 0) that the complete orbit is a piece of the polar axis, penetrating neither 
£=0 nor 77 = 0. 

In the rest of the paper we exclude the possibility of orbits that lie only on the polar axis, 
thus ruling out, e.g., the case of a missile fired vertically from either pole. We then consider 
only those orbits whose completions either reach the plane 77=0 or reach at least one genuine 
extremum 771 for which y = 0, t?^0. 

From (50) it now follows that 

S[^) = ±(e+l)- 1 [al+(k+2^+2a 1 cr)(^+l)]" (53.1) 

^(7 7 ) = ±(l-7 7 2 )- 1 [-ai+(-t+2 M 57 7 +2a 1 cV)(l->7 2 )] i ', (53.2) 

where by (13), (16), and (18) the upper sign is to be used whenever the corresponding co- 
ordinate is increasing and the lower sign whenever it is decreasing. Also 

#(*)=<*. (19) 

From (53) and (19) it then follows that 

^=a 3 0±J^£ 2 +l)- 1 [ai+(t+2 M ^+2a 1 C 2 aa 2 +l)]^ 

±Jj(l-r7 2 )-H-a!+(-*+2^+2a 1 cV)(l-^ 2 )]^+const. (54) 

Before defining the lower limits & and 771, it is well to note that eq (54) is somewhat deceptive 

in appearance, because of the presence of the "small" distance c in terms that are clearly very 

important. Let us now remove this difficulty, as follows. Replace the variable £ by the 

variable p, where 

P^el (55) 

By (41) it then follows that 

p2 =r 2_( 1 _ |? 2 )c 2 (56>1) 

P rj=r sin 6 (56.2) 

¥+1 = (p 2 +c 2 )c- 2 , (56.3) 

so that p«r on the physically realizable part of any complete orbit. On substituting (55) 
into (54) and dropping the constant term, we obtain 

S=a 3 <t>± f * (p 2 +c 2 )- 1 [c 2 ai+(^+2 / xp + 2a lP 2 )(p 2 +c 2 )]irfp 

± f\l-v z )- 1 [-^+(-k+2„8 v +2a 1 cW)(l-V 2 )} i d v , (57) 
where 



ft=<&. (57.1) 

If the complete orbit intersects the plane 77=0, we choose the lower limit 771 to be zero, 
corresponding to 0=0, in accord with the usual canonical treatment [8] of a bound orbit for the 
case V=—fir~ 1 ; this appears to be a good choice for an unbound orbit also. If the complete 
orbit does not intersect the equatorial plane, we choose 77! as an extremum of 77 on the complete 
orbit. It will be a zero of the 77-integrand satisfying 77^ -1. 
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To choose a value for p b we reason as follows. The p-integrand of (57) has the factor 

F( P )=cW+(p 2 +c*)f(p), (58.1) 

where 

/( P ) = 2a lP 2 -f2/xp+£. (58.2) 

The appropriate real zeros of F(p) will be almost equal to those of /(p), which always has two 
in any physically realizable case. If «i=^0, one zero of f(p) will be positive and the other 
negative (approaching— °° as a x approaches zero). If «i<C0, both zeros of/(p) will be positive. 
Any other real zeros of F(p) will occur for inappropriate values of p, of the order c. Thus for 
unbound orbits, we choose pi to be the largest positive zero of F(p), this being the only one 
accessible to complete unbound orbits. For bound orbits, we choose p x to be the next to the 
largest positive zero of F(p), equal approximately to the smaller zero of /(p). In either case 
Pi is then the minimum value of p dynamically accessible to the complete orbit; the correspond- 
ing spheroid may or may not lie within the earth. (Any value of p less than p y as just defined 
is dynamically inaccessible to the complete orbit in the field V(%,r]), even if the planetary matter 
could be removed without change in V(^rj) : to reach such values the orbit would have to 
pass through regions in which F(p) is negative and p$ complex.) Of course p\^>ll and thus 
outside the earth for any ordinary satellite in a bound orbit. 

Whenever the plane r?=0 is intersected by the complete orbit, we have seen that k^ — a*, 
so that we can then introduce a real a 2 such that 



where we can assume 

without loss of generality. Then 



k=-al (59.1) 

a 2 >0 (59.2) 

al^al (59.3) 



If we put 5 = and c = in (57) and then use p = r, 77= sii 1 6, which would (hen follow from (56), 
and also use (59.1), we find that (57) reduces [8] to the correct expression for the action func- 
tion for the case V=— /*r _1 . Thus the terms in (57) that contain c 2 or <5 give rise to deviations 
from a conic-sectional orbit. 

From (13.3) and (19) it follows that the constant a z in (57) is the ^-component of angular 
momentum. The constant k has no name, but is related to the magnitude of the total angular 
momentum, which is, however, not conserved in fche non-centra] field produced by oblateness. 
To see the connection, note that whenever k^—a 3 2 , we can write (59.1) and (59.3). Then 
when 8 and c 2 are both allowed to approach zero in (57), a 2 approaches the magnitude of the 
total angular momentum, which is then conserved. (This result follows most easily from the 
fact already mentioned, viz, that (57) then reduces to the correct action function for the case 
V=-»r~\) 

9. Kinetic Equations of Motion 
Implicit equations for p, 77, and <f> as functions of t are now given by 

ZS/da^t+fa (60.1) 

dS/dk=p' 2 (60.2) 

bS/da 3 =p 3 (60.3) 

where t is the time and the jS's are constants. Since the lower limit p x in (57) is a function of 
a u k, and <* 3 , the derivative of the p-integral with respect to any of them is the sum of an in- 
tegral and a term like TlpCp^dp^dai. Here I p (p) is the p-integrand, which vanishes for pi, 
so that this second term vanishes. If the complete orbit intersects the plane 77=0, so that 771 
can be taken to be zero, on second term arises in a derivative of the 77-integral; if not, rji is to be 
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taken as an extremum value of 77 on the complete orbit and the second term again vanishes, 
since rji is then a zero of the 77-integrand. 
The kinetic equations are thus 



± fV[^!+(&+2 M p+2a lP W+^ 

J p x Jv 1 

(61.1) 

±i( P {c*al+(k+2,xp + 2a 1 p 2 )(p i +C i )]-id P T$ f''[-«!+(-^+2^+2a I cV)(l-'/ 2 )]- i ^ = ^, 

(61.2) 
tf-ic^J V+c 2 )~ V«!+ (k+2„p + 2a lP 2 ) (p 2 +c 2 )]-idp 

T« 3 [\l-^- 1 [-al+{-k+2^+2a l cW){l-^]-'dr ] = ^. (61.3) 

If we restrict considerations to those cases where the plane 17 = is intersected hy the 

complete orbit, then k^a\. We can then substitute k= — a 2 2 ^—a 2 3 in (57) and replace 
(60.2) by 

dS7c)a 2 = |3 2 . (62) 

From (57), (60.1), (60.3), and (62) there then follow the kinetic equations 

±£ p 2 [c 2 al+{-a\+2 t ip+2a l p 2 ){p 2 +c 2 )]-Up±c 2 ^ V 2 [-al+(al+2„S V 

+ 2a i cy)(l- V 2 )]-idr, = t + p 1 , (63.1) 

T« 2 P [c 2 al+(-«!+2 ti p+2a 1 p 2 )(p 2 +c 2 )]-idp±a 2 i" [-al+(al+2 l xS V +2a 1 c 2 r, 2 )(l-v 2 )]-klr,=p 2> 
J i>i Jo 

(63.2) 

<P±C 2 a 3 f" (p 2 +C 2 )- 1 [c 2 al+(-«I+2MP+2« I p 2 )(p 2 +c 2 )]-Wp 

J Pi 

T a^\l-^)~'[-al+{al+2 l ,by ] +2a l cW)il- (63.3) 

Future papers will be devoted to evaluating the integrals in these equations, solving them 
simultaneously for p, 77, and <j> as functions of time, and deriving the astronomical results that 
they imply, both for bound orbits and unbound orbits. 

10. References 

[1] J. P. Vinti, Phys. Rev. Letters 3, No. 1 (July 1959). 

[2] R. H. Merson and D. G. King-Hele, Nature 182, 640 (Sept. 1958). 

[3] H. Jeffreys, The earth, 3d ed (Cambridge Univ. Press, Cambridge, England, 1952). See p. 127 at bottom; 

p. 129, eq (5); and p. 132, eq (13). 
[4] D. G. King-Hele and R. H. Merson, Nature 183, 881 (March 1959). 
[5] J. A. O'Keefe, A. Eckels, and R. K. Squires, Science 129, 565 (Feb. 1959). 
[6] J. A. O'Keefe, Astron. J. 62, 265 (1957). 
[7] P. M. Morse and H. Feshbach, Methods of theoretical physics, vols. I and II (McGraw-Hill Book Co., 

Inc., New York, N.Y., 1953). 
[8] W. M. Smart, Celestial mechanics, p. 143 (Longmans, Green & Co., New York, N.Y., 1953). 

Washington, D.C. (Paper 63B2-12) 



116 



